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ABSTRACT 

We explore the effects of dissipationless merging on the Fundamental Plane of 
elliptical galaxies by using a N-body code based on a new, high performance numeri- 
cal scheme (Dehnen 2002). We investigate the two extreme cases of galaxy growth by 
equal mass merging and accretion of small stellar systems; in a subset of simulations 
we also consider the presence of dark matter halos around the merging galaxies. Cu- 
riously, we found that the Fundamental Plane is preserved by major merging, while 
in the accretion scenario its edge-on thickness is only marginally reproduced, with 
substantial thickening in the case of merging with low angular momentum. We also 
found that both the Faber- Jackson and Kormcndy relations are not reproduced by 
the simulations, in accordance with the results of a preliminary analysis based on a 
simple application of the virial theorem. Finally, we discuss the implications of our 
results for the origin of the Mbh-co and Magorrian relations. We found that dissipa- 
tionless merging is unable to reproduce the A/bh-co relation, if the black hole masses 
add linearly (while the Magorrian relation is nicely reproduced); on the contrary a 
black hole merging with substantial emission of gravitational waves reproduces the 
Mbh-0o relation but fails at reproducing the Magorrian relation. We argue that our 
results strongly point towards a major role of dissipation in the formation of early-type 
galaxies and in the growth of their central supermassive black holes, thus supporting 
the idea of a link between galaxy formation and QSO activity. 
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1 INTRODUCTION 

Two distinct scenarios have been proposed to describe the 
formation of elliptical galaxies (Es): roughly speaking, ac- 
cording to the monolithic scenario Es form at high red- 
shift in dissipative collapses (see, e.g., Eggen, Lynden-Bell & 
Sandage 1962, Larson 1975), while in the hierarchical merg- 
ing scenario spheroidal systems are the end-products of sev- 
eral merging processes of smaller galaxies, the last major 
merger taking place in relatively recent times, i.e. at z<l 
(see, e.g., White & Rees 1978, Kauffmann 1996, Cole et 
al. 2000). In particular, this second picture is supported by 
some observational data suggesting that a fraction of red 
galaxies in clusters at intermediate redshift are undergoing 
merging processes: these galaxies could be the progenitors of 
present-day early-type galaxies (van Dokkum et al. 1999). 

In any case, it is well known that Es satisfy many scaling 



relations (more or less tight) that must be accounted for by 
any proposed formation scenario: for example the Faber- 
Jackson relation (hereafter FJ, Faber & Jackson 1976), the 
Kormendy relation (Kormendy 1977), the Mg 2 -cro relation 
(see, e.g., Burstein et al. 1988, Bender, Burstein & Faber 
1993) the color-magnitude relation (Bower, Lucey & Ellis 
1992), and the very tight Fundamental Plane (hereafter FP; 
Djorgovski & Davis 1987, Dressier et al. 1987). In addition, it 
is now widely accepted that central supermassive black holes 
(hereafter BHs) are a common characteristic of spheroidal 
stellar systems (see, e.g., Kormendy & Richstone 1995, van 
der Marel 1999, de Zeeuw 2001); it has also been found that 
their mass (Mbh) is linearly proportional to the mass of 
the host spheroid (Magorrian et al. 1998), and it is even 
more correlated with the spheroid central stellar velocity 
dispersion (Gebhardt et al. 2000, Ferrarese & Merritt 2000). 
Due to their relevance in this work, here we briefly review 
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the main properties of the FP, FJ, Kormendy and Mbh-o"o 
relations. 

The FP of Es is a scaling relation between three of their 
basic observational properties, namely the circularized effec- 
tive radius (R) e = V o c b c (where o e and b e are the major and 
minor semi-axis of the effective isophotal ellipse), the cen- 
tral velocity dispersion ao, and the mean effective surface 
brightness (I) e =LB/2n(R) c 2 (where Lb is the luminosity of 
the galaxy, for example in the Johnson B-band) . A parame- 
terization of the FP, that we adopt in this paper, has been 
introduced by Bender, Burstein & Faber (1992, hereafter 
BBF), by defining the three variables 



k 2 = 



logcrg + log(fl)c 

logqg + 21og(I)c-log(-B)c 

loggg-log(/)o-lQg<fl)o 

v/3 



(1) 
(2) 
(3) 



In particular, when projected on the (ki,ks) plane, the FP 
is seen almost edge-on and it is considerably thin, while the 
distribution of galaxies in the (fei, fe) plane is broader. For 
example, Virgo ellipticals studied by BBF are distributed on 
the (k\,kz) plane according to the best-fit relation 



0.15fci +0.36, 



(4) 



(when adopting respectively, kpc, km s _1 and Lb© pc~ 2 as 
length, velocity and surface brightness units), with a very 
small l-cr dispersion of rms(£;3) ~ 0.05 over the whole range 
spanned by the data, 2.6 < ki < 4.6 (0.75 < h A < 1.05). This 
translates in a scatter of ~ 15% in {R) c for fixed Lb and ao- 
The slope and the scatter of equation (4) are usually called 
"tilt" and "thickness" of the FP, respectively. By combining 
equations (1) and (3) with equation (4) the FP relation is 
obtained in terms of the observables, and is in good agree- 
ment with the FP derived (in the same photometric band 
and for a much larger galaxy sample) by j0rgensen, Franx 
& Kjaergaard (1996)*. 

Strictly related to the FP are the less tight FJ and Kor- 
mendy relations. These scaling relations have been often in- 
terpreted just as projections of the FP with no additional in- 
formation, and, for this reason, their constraining power on 
galaxy formation scenarios has been underestimated. How- 
ever this is misleading, because these two relations describe 
how Es are distributed on the FP and so, even if character- 
ized by a larger scatter than the edge-on view of the FP, 
they contain important information on the galactic prop- 
erties. The proposed original form of the FJ relation was 
Lb oc ctq , with n « 4, while Davies et al. (1983) found that 
the double-slope fit 

life * °- 23 (scK^f 4+0 - 62 (m£^T 2 ^ 

provides a better description of the data. Note that for small 
galaxies and bulges (ao < 170 kms -1 ) the exponent is ~ 2.4, 



* As it is well known, the FP tilt depends on the adopted pho- 
tometric band. Curiously, when using the fc-space, the coefficient 
in the K-band is only slightly different from 0.15, with a reported 
value of 0.147 (Pahre, Djorgovski & dc Carvalho 1998b). 



considerably smaller than 4; this situation is also reflected 
by the fit of Dressier et al. (1987) to the galaxies of Virgo 
cluster, Lb oc o"o' 5 - Due to indeterminacy of the exact value 
of n, in the rest of the paper we use 3.5 and 4 as represen- 
tative values of the FJ exponent. 

The total luminosity of spheroidal systems also corre- 
lates with their length scale as measured by (R) e : in fact, 
the Kormendy relation can be written in the form 



(R) e oc Lq, 



(6) 



where the exponent a is strongly dependent on the galaxy 
sample used, and is found in the range 0.88 <a< 1.62 
(Ziegler et al. 1999). The latest estimates seem to con- 
verge to a value a ~ 0.7 or less, as a function of waveband 
(Bernardi et al. 2003). 

In the last years, another important scaling relation has 
been added to the list, the so-called Mbh-oo relation (see, 
e.g., Gebhardt et al. 2000; Ferrarese & Merritt 2000, Merritt 
& Ferrarese 2001, Tremaine et al. 2002): 



M B H OC <Jq, 



(7) 



where ao is the projected central velocity dispersion of the 
parent galaxies (or within (R} c ), and the exact value of a 
is still matter of debate, ranging between 4 and 5. An im- 
portant characteristic of this relation is its extremely small 
scatter, consistent with measurements errors only, so that 
equation (7) can be considered a "perfect" relation. 

In the context of galaxy formation studies, a natural 
question to ask is how well dissipationless merging is able 
to produce and maintain these relations. We focus on dis- 
sipationless merging since it is simpler to be modeled with 
respect to gas rich merging, and because in this way we can 
check whether the contribution of gas dissipation is required 
in the merging scenario of Es formation. In addition, the 
analysis of dissipationless merging should indicate whether 
the merging between gas poor systems observed at z < 1 
(van Dokkum et al. 1999) is a common phenomenon in Es 
lifetime, or a rare event. Among others, existence of serious 
problems encountered by the picture of Es formation mainly 
driven by dissipationless merging to reproduce their scaling 
relations was pointed out by Ciotti & van Albada (2001; 
hereafter CvA). In their phenomenological approach, just by 
combining the observed (edge-on) FP and the Mbh-co re- 
lations, they demonstrated that in a dissipationless merging 
scenario, and under reasonable assumptions for the addition 
of BH masses during the merging, elliptical galaxies forced 
to he on the FP and to satisfy the Mbh-co relation would 
have effective radii exceedingly larger than those observed 
in real Es. 

Here we follow a complementary approach, and the 
main goal is verifying, by using high resolution N-body nu- 
merical simulations of one and two-component galaxy mod- 
els, whether the end-products of merging of galaxies, ini- 
tially lying on the FP, lie on the FP (as well as on the other 
scaling relations) or they fail in some other respect (see also 
Pentericci, Ciotti & Renzini 1995; Capelato, de Carvalho 
& Carlberg 1995; Bekki 1998; Evstigneeva, Reshetnikov & 
Sotnikova 2002; Nipoti, Londrillo & Ciotti 2003a). 

Numerical simulations are needed since, owing to possi- 
ble structural and/or dynamical non-homology and to pro- 
jection effects, one cannot predict with simple theoretical 
arguments how ao and (R) e evolve as consequence of merg- 
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ing: however, rough indications on their behavior can be 
derived from the study of the evolution of the virial velocity 
dispersion ay and the virial radius ry^ . In fact, as a con- 
sequence of the virial theorem and the conservation of the 
total energy, in the merging of two galaxies with masses M\ 
and M2 and virial velocity dispersions crv,i and ffv,2, the 
virial velocity dispersion of the resulting galaxy (in case of 
no mass loss and negligible kinetic and interaction energies 
of the galaxy pair when compared to their internal energies) 
is given by 

2 _ M lf r v>1 + M 2 a v>2 
av < 1+2 ~ M 1+ M 2 • (8) 

It follows that ov, i+2< max(crv,i, ov, 2), i.e., the virial ve- 
locity dispersion cannot increase in a merging process of the 
kind described above. Under the same hypotheses, the virial 
radius rv,i+2 of the resulting galaxy is given by 

(Mi + M 2 f _ Ml | Ml 
?*v,i+2 rv,i r v ,2 ' 

where rv,i an d ry,2 are the virial radii of the progenitors. 
Identity (9) implies that 

rv,i+2> min(rv,i,fv,2), i.e., the virial radius cannot de- 
crease in a merging process of the kind described above. 

From equations (8) and (9) it follows that, in the highly 
idealized scenario of a merging hierarchy based on identical, 
one-component seed galaxies characterized by ov,o, fv,o and 
Mo, we expect ay = <jv,o and rv = (M/Mo)rv,o, indepen- 
dently of the merging sequence. Thus, by qualitatively as- 
suming that o~o ~ ay and (R} c ~ rv, one should conclude 
that the FJ and Kormendy relations are not consistent with 
the scenario depicted above. Now, it is well known that for a 
large variety of mass models there is a good correlation be- 
tween tm (the half-mass radius of the galaxy, strictly related 
to (R) e ) and ry, with the proportionality constant showing 
little dependence on the particular density profile^ (see, e.g., 
Spitzer 1969, Ciotti 1991), and so, even in presence of struc- 
tural non-homology in the merging end-products, the "pre- 
diction" of a linear growth of (R) e with M should be quite 
robust. A more problematic situation (that can be properly 
addressed only with numerical simulations) arises when con- 
sidering the relation between cry and ao: at variance with 
rm/ry, structural non-homology can strongly affect the ra- 
tio go/ cry (see, e.g., Ciotti, Lanzoni & Renzini 1996; Bertin, 
Ciotti & Del Principe 2002). In addition, a systematic vari- 
ation of orbital anisotropy with mass (the so-called dynam- 
ical non-homology) can produce an increasing o~o even at 
constant ay (although Nipoti, Londrillo & Ciotti 2002, here- 
after NLC02, showed that the whole tilt of the FP cannot 
be ascribed to orbital anisotropy effects only). 

Here, with the aid of high-resolution N-body simula- 
tions of one and two-component galaxy models, we investi- 
gate the effect of subsequent generations of merging on the 
galaxy scaling relations. Our study explore two extreme sit- 
uations, namely the case of major merging, in which equal 



t By definition in a one-component galaxy cry = 2T/M and 
ry = —GM 2 /U, where T and U are the total kinetic and the 
gravitational energy of the galaxy, respectively, 
t For example in truncated power-law spherical models ry /r m = 
(5 - 2 7 )/(3 - 7) and so 2 < r v /r M <2.1 for < 7 < 2. 



mass galaxies are involved at each step of the hierarchy, and 
the case of accretion, in which a massive galaxy increases 
its mass by incorporating smaller galaxies: the evolution of 
more realistic, "mixed" merging histories should be brack- 
eted by our simulations. This paper is organized as follows. 
In Section 2 we describe the properties of the galaxy mod- 
els and the adopted numerical methods. In Section 3 we 
present the results, while in Section 4 we specifically discuss 
the consequences of the findings on the implications for the 
Mbh-to relation. Finally, in Section 5 the main conclusions 
are summarized. 



2 NUMERICAL METHODS 
2.1 Galaxy models 

As initial conditions for the first generation of merging we 
use spherically symmetric one and two-component Hern- 
quist density distributions (Hernquist 1990; Ciotti 1996, 
1999). This choice is motivated by the fact that the Hern- 
quist model is both a reasonable approximation (when pro- 
jected) of the i? 1//4 law (de Vaucouleurs 1948), and it is also 
similar to the Navarro, Frenk & White (1996) profile, thus 
providing a sufficiently realistic description of the galaxy 
stellar and dark matter (DM) density distributions. The 
density, mass, and (relative) potential profiles of the stel- 
lar component of the Hernquist model are given by 

M,r c 



P * {r) ~ 27rr(rc + r) 3 ' 

M«(r) _ / r \ 2 
M, ~~ \r c + r) ' 
GM, 



(10) 



(11) 



(12) 



r c + r 

where M» is the total stellar mass. In the two-component 
models, the DM halo is described by p b , M h and * n profiles 
of the same family of equations (10)-(12), where now r n = 
0r c and Mh = /J,M*. At variance with the models used in 
NLC02, here we restrict to globally isotropic models: the 
distribution function (DF) of the stellar component is then 
given by 



d 



I 



dp* d* T 



y^TT 2 d£ J d* T yJE - 



(13) 



where £ — $t — v 2 /2 is the relative (positive) energy, v is 
the modulus of the velocity vector, and $t = + *h is 
the relative total potential; an analogous identity holds for 
the DM halo, where ph substitutes p«. 

According to the given definitions, the one-component 
models are completely determined by the two physical scales 
M* and r c , and in the numerical simulations M* , r c and Td yn 
are adopted as mass, length and time scales. In particular 
Tdyn is the half-mass dynamical time, defined as 



'dyn 



3tt 



8.33 



GM, 



(14) 



16Gp M 

where p M = 3M* /87rrM is the mean density inside the half 
mass radius vm- Finally, the velocity scale is given by 



' dyn 



0.12 



GM, 



(15) 
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The two-component Hernquist models are characterized by 
four quantities, the two physical scales M» and r c , and the 
two dimensionless parameters fi and (5. In order to transform 
the results of the numerical simulations in physical units we 
express M* in 10 10 Mq, r c in kpc, Lb in 1O 1O Z/B0, and the 
constant stellar mass-to-light ratio T* = M»/Lb in solar 
units. For a detailed description of the numerical realization 
of the initial conditions, see Section 2.1.3 of NLC02. 

2.2 Initial conditions 

2.2.1 Equal mass merging 

The first generation of the equal mass merging hierarchy is 
obtained by merging a pair of identical, spherically symmet- 
ric and isotropic Hernquist models (the "zeroth order" seed 
galaxies), while the successive generations are obtained by 
merging pairs of identical systems obtained by duplicating 
the end-product of the previous step. We followed the evo- 
lution of the hierarchy of 5 steps of head-on mergers (that 
would correspond, in absence of mass escape, to a mass in- 
crease of a factor 32); we also explored 3 steps of the head-on 
merging hierarchy whose seed galaxies are two-component 
Hernquist models (in which we assume a DM halo more 
massive, fj, = 3, and less concentrated, (3 = 2, than the 
stellar component), and, finally, 3 steps of the hierarchy of 
encounters of one-component models with non zero (orbital) 
angular momentum. 

In the assignment of the initial conditions for the head- 
on encounters the two galaxies are initially placed on an 
orbit characterized by vanishing relative energy and angular 
momentum, i.e., they would have a null relative velocity 
while at infinite relative distance. In practice, at the time 
t = 0, we place the two galaxies at a distance d rc i ~ 3(r)9o 
(where {r)go is the angle-averaged radius enclosing the 90% 
of the total mass M of each galaxy), and, neglecting the 
effects due to tidal forces, we assign to th eir center s of mass a 
relative velocity with modulus t> rc i = 2y/ GM/d le \. The end- 
products of each merging are non-spherical, and so, when 
exploring a successive merging, their mutual orientation at 
the beginning of the new simulation is randomly assigned. 
In the simulations with orbital angular momentum, we still 
consider vanishing relative energy and d rc i defined as above, 
and we assume as impact parameter the sum of the angle- 
averaged half-mass radii of the merging galaxies. 

2.2.2 Accretion 

In the case of the merging hierarchy in which the test galaxy 
grows by accretion of smaller systems {accretion), the seed 
galaxy is a one-component Hernquist model, and the first 
merging event is identical to that in the equal mass merging 
case. In the second step, however, the first end-product (of 
mass ~ 2M,) merges again with a seed galaxy, and so on. As 
a consequence, the mass ratio between the infalling stellar 
system and the test galaxy decreases approximately as 1/n, 
where n is the step in the accretion hierarchy. Due to the 
higher computational cost required by accretion simulations 
with respect to equal mass merging in order to reach the 
same mass increase, in the former case we limit to 15 steps 
(of head-on merging) and to 9 steps (of merging with angu- 
lar momentum), for a putative mass increase of a factor 16 



and 10, respectively. As in the case of equal mass merging 
simulations, the relative initial positions and velocities of the 
infalling satellite and of the test galaxy are chosen so that 
at t = they are on parabolic orbit, and the mutual orien- 
tation of the two galaxies is assigned randomly. In practice, 
as initial relative distance we adopt d rc i ~ (r)go,i + 2(r)go,2, 
where the subscript 2 refers to the smaller galaxy. Again, in 
case of non zero angular momentum, the impact parameter 
is given by the sum of the angle-averaged half-mass radii. 

2.3 The numerical code 

For the simulations we used a new, fast and accurate, 
parallel N-body code. Recently, a numerical scheme has 
been introduced to simulate collisionless N-body systems 
by Dehnen (2000, 2002), with a substantial improvement 
in force computation over standard tree-based codes. By 
combining in an original way the Fast Multipole Algorithm 
(see, e.g., Greengard & Rokhlin 1987, 1997) with the Barnes 
& Hut (1986) tree data structure, Dehnen scheme achieves 
exact momentum conservation with an effective O(N) oper- 
ational complexity (for particle numbers AT>10 ). 

By taking advantage of these improvements, we have 
implemented the Dehnen fast force solver (FalcON) in a 
parallel Fortran-90 N-body code named FVFPS (Fortran 
Version of a Fast Poisson Solver; Londrillo, Nipoti & Ciotti 
2003). The FVFPS code is based on domain decomposition, 
and uses the MPI routines for data communication. A mass 
dependent opening parameter 6 — 9(M), that assures faster 
performances and a more uniform error distribution, is also 
employed (Dehnen 2002). Time integration is performed by 
a standard leap-frog algorithm, with (uniform) time step At 
determined adaptively. Thanks to effective O(N) scaling of 
the force computations, the FVFPS code assures a speed-up 
of a factor ~ 10 over existing codes, in particular over the 
Springel, Yoshida & White (2001) GADGET code. In the 
present version, the FVFPS code is characterized by three 
preassigned parameters: 

(i) The minimum value of the opening parameter # m i n , 
associated with the total mass. We adopt # m in = 0.5, result- 
ing in mean relative force errors of ~ 0.3%, and 99 percentile 
errors of ~ 3%. 

(ii) The softening parameter e (i.e., the softening length 
expressed in units of r c ), dependent, according to literature 
indications (e.g., Merritt 1996, Athanassoula et al. 2000, 
Dehnen 2001), on N and on the specific density distribu- 
tion profile. 

(Hi) For the initial time step we assume At ~ T dyn /100. 
The time step is allowed to vary adaptively as a function of 
the maximum particle density, but is kept the same for all 
the particles. 

For the stellar distribution of the seed galaxies we use 
as a rule AT* = 16378 particles for each galaxy, and in 
the two-component cases the DM distribution is made of 
Num = 49152 particles: with this choice, halo and stellar 
particles have the same mass. As a consequence of the merg- 
ing hierarchy, the number of particles in the simulations in- 
creases with the galaxy mass: in the final merging of (one 
and two-component) equal mass systems and of the accre- 
tion scenario, the total number of particles involved is of the 
order of 5.2 x 10 5 and 2.6 x 10 5 , respectively. 

We followed the dynamical evolution of each merging 
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-l o 
log r/r M 

Figure 1. Anglc-avcragcd density as a function of radius (nor- 
malized to the half-mass radius tm) of the end-product of the 
head-on merging of two one-component Hernquist models, ob- 
tained with GADGET (empty triangles) and with the FVFPS 
code (solid circles). 
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Figure 2. Scrsic best fit parameter m vs. total stellar mass of 
the end-products at stage i of the merging hierarchy. Equal mass 
mergers arc shown as solid triangles and squares (one-component 
galaxies), and stars (two-component galaxies); empty triangles 
and circles represent the accretion hierarchies. Triangles corre- 
spond to simulations with non zero orbital angular momentum. 



event up to the virialization of the resulting system, which 
is usually reached on a time scale shorter than 50Td yn after 
the first encounter between the two galaxies: for example 
the virialization time for the first step in the one-component 
(head-on) equal mass merging hierarchy is of the order of 
2Gyr having assumed M» = 1O 1O M0 and r c = lkpc. In 
the following we call "end-product" of a simulation the sys- 
tem made of bound particles after the virialization. For the 
determination of the intrinsic and "observational" proper- 
ties of the end-products, we followed the same procedure 
described in NLC02. In particular we measured the end- 
product intrinsic axis ratios c/a and b/a (where a, b and 
c are the major, intermediate and minor axis of the asso- 
ciated inertia ellipsoid), their virial velocity dispersion cry 
and angle-averaged half-mass radius rm, and, for several 
projection angles, the end-product circularized effective ra- 
dius (R) e , the central velocity dispersion oo, and the mean 
effective surface brightness (7) e (obtained from the stellar 
mass profile by assuming a constant stellar mass-to-light 
ratio T»). In particular, the central velocity dispersion ao is 
obtained by averaging the projected velocity dispersion over 
the circularized surface brightness profile within a radius of 
{R) e /8- As a general rule we find that discreteness effects on 
the derived values of (R) e and er do not exceed 1 -j- 2% at 
each projection angle. 

We checked the reliability of the numerical simula- 
tions by running a few merging events both with the par- 
allel version of GADGET (with a = 0.02, Ai min = 0, 
Af max = Td y n/100, atoi = 0.05, and e = 0.05 ;see Springel 
et al. 2001 for details) and with our FVFPS code (with 
6» min = 0.5, e = 0.05 and At = r dyn /100). The results 
are in remarkable agreement: the end-products in three- 
dimensional shape, angle-averaged half-mass radius, den- 
sity and velocity dispersion profiles are practically indistin- 



guishable. With both codes the total energy is conserved 
within 1% over 100 dynamical times. As an illustrative ex- 
ample, in Fig. 1 we show the density profiles of the end- 
products of the head-on merging of two Hernquist mod- 
els (GADGET, empty triangles; FVFPS, solid circles); note 
how "minor" details, as the slope change in the density pro- 
file at r ~ 10rM, are nicely reproduced. All the simulations 
were run on a Cray T3E, and on a IBM Linux Cluster (with 
a number of processors from 4 up to 32, depending on the 
number of particles of the simulation). 



3 THE RESULTS 

3.1 Equal mass merging 

Before presenting the effects of equal mass merging on the 
FP and on the other scaling relations, we briefly describe the 
induced modifications of the internal structure and dynam- 
ics of one and two-component galaxy models. In general, 
the end-products are triaxial systems with axis ratios in the 
range 0.5 < c/a < 0.7 and 0.7<6/a<0.8 (where a, b and c 
are the major, intermediate and minor axis, respectively), 
in accordance with observed ellipticities of Es. For compar- 
ison with real galaxies, we also fitted (over the radial range 
0.1 < R/{R) e < 10) the projected stellar mass density profile 
of the end-products with the Sersic (1968) R 1/m law: 



I(R) = Io exp 



(I 



l/m 



(16) 



where b(m) ~ 2m-l/3+4/405m (Ciotti & Bertin 1999), and 
in Fig. 2 we plot the best fit parameter m as a function of the 
total mass of the systems, for one-component (solid squares 
and triangles) and two-component (stars) models. Clearly, 
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the fitted quantities m and {R) e depend on the relative ori- 
entation of the line-of-sight and of the end-products of the 
simulations: the two points for each value of the mass in Fig. 
2 show the range of values spanned by m when projecting 
the final states along the shortest and longest axis of their 
inertia ellipsoids. We note that higher values of the best fit 
parameter m correspond to more massive systems, a trend 

similar to that observed in Es§ (see, e.g., Caon, Capaccioli & 
D'Onofrio 1993, Prugniel & Simien 1997, Bertin et al. 2002). 
Note also that the values of m are in the same range obtained 
from observations, in fact we found 2 < m < 11. (However, as 
pointed out to us by the Referee, from Fig. 2 one could ar- 
gue that m ~ 4 galaxies should have experienced at most 
one major merger event in their life.) 

Following the discussion in the Introduction, in Fig. 3 
(top panel) we show the relation between the virial veloc- 
ity dispersion and the total stellar mass of the mergers at 
each step of the merging hierarchy. In one-component mod- 
els (solid squares and triangles) a modest increase of a\r 
with M» is apparent: for the last model in the hierarchy the 
mass is increased by a factor ~ 24.5 (instead of the maxi- 
mum possible value of 32), while the ov is increased by a 
factor ~ 1.16 only. This small increase of try with respect 
to the expectation of equation (8) (i.e., ov = const) can be 
explained with a simple generalization of the highly ideal- 
ized situation formalized there. In fact, if during a merging 
between two galaxies with mass Mi , M 2 and virial velocity 
dispersion <7v,i, <rv,2 a mass AM is lost with mean velocity 
w c j, then the virial velocity dispersion of the end-product 
ov,i+2 is expected to be 



2 _ Mio- Vil + M 2 o- V;2 + AM<j 

°V,l+2 — ~~ ~~ ~~ ~ 



Mi + M 2 - AM 



(17) 



in particular <rv increases as a consequence of the mass (and 
associated kinetic energy) loss. We found that in the sim- 
ulations the mass lost during each merging never exceeds 
AM/(Mi +M2) ~ 0.06 and, consistently, the corresponding 
increase of av is 2 4- 5% (Fig. 3, top panel). 

In Fig. 3 (bottom panel) we plot the half-mass radius 
as a function of the stellar mass of the mergers. As for ov, 
also for this quantity the simple virial expectation of a linear 
growth of rv (and of ru) with M (dotted line) is apparent. 
The "jump" from the initial condition position to the first 
merger is due to the significant change in the galaxy density 
structure (as revealed by the change of values of m in Fig. 
2); the successive generations of merging are characterized 
by more similar density profiles, and, correspondingly, the 
points in Fig. 3 move parallel to the dotted line. 

How does the presence of massive DM halos affect these 
results? As anticipated, in this case we investigated 3 steps 
of the (head-on) equal mass merging hierarchy by using 
two-component galaxy models. Due to its observational rel- 
evance, we focus here on the description of the properties of 
the stellar component only. As in the one-component case, 
the end-products are triaxial, and with the axis ratios in 
the same range. The best fit Sersic parameter m still in- 



§ Curiously, NLC02 found that in the case of end-products of 
unstable galaxy models m decreases for increasing radial orbital 
anisotropy in the initial conditions. Note also that in NLC02 and 
Londrillo ct al. 2003 we used a radial range 0.1 < -R/(R) e <4. 



creases with mass (stars in Fig. 2). In the two-component 
simulations, owing to the presence of the DM halo, there are 
not arguments as simple as those used in equations (8), (9) 
and (17) in order to predict the effects of merging on ov 
(defined as 2T*/M„) and ryi of the stellar component. How- 
ever, it is apparent from Figs. 2 and 3 that the presence of a 
DM halo or angular momentum in the initial conditions does 
not modify the general trend of m, ov and rm at increasing 
mass. 



3.1.1 Fundamental Plane 

In Fig. 4 we plot the results of equal mass merging sim- 
ulations in the principal planes (ki,k:j) and (fci,fc 2 ), where 
one-component galaxy models resulting from head-on merg- 
ers and from merging with angular momentum are identified 
by solid squares and triangles, respectively; two-component 
models are represented by stars. The progenitor of the merg- 
ing hierarchy (the point without bar) is placed on the edge- 
on FP by choosing its luminosity Lb — 3 x 10 9 Lbq and 
stellar mass-to-light ratio T* = 5, so that equation (4), 
represented by the solid line, is satisfied. Consistently with 
the adopted dissipationless scenario, the value of T» is kept 
constant during the whole merging hierarchy. Due to the loss 
of spherical symmetry of the end-products of the merging 
simulations, their coordinates depend on the line-of-sight 
direction; however, being the luminosity (mass) of each end- 
product fixed, variations of fci, k'2, hi due to projection ef- 
fects are not independent. In fact, the fci and k$ coordinate 
of a galaxy of given luminosity are linearly dependent as 



fc 3 



(18) 



this dependence is reflected by the straight lines in Fig. 4 
(top panel), each of them representing the range spanned 
in the (fci, fcs) space by a given end-product when observed 
over the solid angle. As in NLC02, in order to quantify the 
deviation of a model from the FP, we use the vertical dis- 
tance Sk$ = \k3 — 0.15fci — 0.36j from the point to the FP 
itself. It is interesting to note that for all the end-products 
the projection effects translate into a 5ks, of the same order 
of magnitude of the observed FP dispersion. In addition it 
is also apparent from Fig. 4 how £3 and fci increase with 
the galaxy mass consistently with the observed FP tilt and 
thickness. Being T* fixed in our simulations, from Fig. 2 
we conclude that one—component equal mass dissipationless 
merging is able to reproduce (basically by structural non- 
homology) the edge-on FP of elliptical galaxies (see also 
Capelato et al. 1995). 

We also explored the behavior of the models in the 
(fci,fc2) plane (Fig. 4, bottom panel), which represents the 
face-on view of the FP; the region populated by real galaxies 
in this plane is identified by dashed lines. By fixing the scale 
length, we place the first progenitor (with {R) c — 0.7 kpc) 
in a zone of the (fci, fc2) plane populated by low-luminosity 
ellipticals (fci ~ 3, fc 2 ~ 4.5), and for each end-product we 
plot the (fci,fcs) positions for a few random projection an- 
gles. Note that, at variance with the coordinates (fci, £3), 
fci and fc2 are related (for a given end-product) by an ex- 
pression containing (R) c besides Lb, and for this reason in 
the (fci, fc 2 ) plane the effect of projection is to distribute the 
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Figure 3. Top panel: virial velocity dispersion of the stellar component at stage i of the merging hierarchy vs. the total stellar mass 
of the merger. Bottom panel: angle— averaged half-mass radius vs. total stellar mass. Symbols are the same as in Fig. 2. The dotted 
line indicates the relation rjyj oc M; M,^, o an d °V are the stellar mass, the virial velocity dispersion and the half-mass radius of 
the seed galaxy. In the two-component cases, <Ty and rn refer to the stellar component only. Note the different range spanned in the 
ordinate axes in the two panels. 



end-products on two-dimensional regions. One-component 
equal mass mergers (solid squares and triangles) are moved 
by merging towards the bottom right of the (fcijfe) plane, 
roughly parallel to the line defining the zone of avoidance, 
in accordance with the prediction of BBF. The produced 
displacements are very large: it is remarkable that the end- 
products of the last step of the hierarchy are found in a 
position only marginally consistent with the populated re- 
gion in the (fcijfe) space. It is also interesting to note that 
the presence of substantial angular momentum in the ini- 
tial conditions does not change significantly the properties 
of equal mass mergers, both in the (ki,ks) and in the 
plane. 

The results of equal mass merging of two-component 
galaxies are represented in Fig. 4 with stars. Due to the 
presence of DM, in order to place the first progenitor (Lb = 
3 x 10 9 Lbq) on the FP, we assumed a stellar mass-to-light 
ratio T» = 2.3. It is apparent that both in the (fci,/ca) 
and (k\,k2) planes, the behavior of two-component mod- 
els docs not significantly differ from the corresponding one- 
component models. In other words, dissipationless, equal 
mass merging of one and two-component models seems to 
be consistent with the existence of the FP of Es. 



3.1.2 Faber- Jackson and Kormendy relations 

As we have shown in the Section above, dissipationless merg- 
ing of one and two-component galaxy models is surprisingly 
consistent with the existence of the FP, especially consider- 
ing its small thickness when seen edge-on. However, Es do 
follow additional scaling relations as the FJ and Kormendy 
relations, and so here we compare the results of our simula- 
tions with these scaling laws. 

The solid line in the upper panel of Fig. 5 represents the 
FJ relation in the form Lb oc a\ together with its scatter 
((51oguo — 0.1, Davies et al. 1983, dashed lines); the dot- 
ted line represents the steeper fit (Lb oc cxq ) derived by 
Dressier et al. (1987). The end-products of (head-on) equal 
mass merging of one-component galaxies (solid squares) 
have a <jo lower than that predicted by the FJ relation for 
the given mass increase, and the discrepancy is stronger in 
case of merging with angular momentum (solid triangles). It 
is interesting to compare this plot with the upper panel of 
Fig. 3: it is then clear how the failure of FJ can be directly 
traced to the constancy of cry and to the significant dynami- 
cal homology of the mergers. This latter point again suggest 
that structural non-homology, more than dynamical non- 
homology is at the origin of the FP in the explored merging 
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Figure 4. Top panel: the merging end— products in the (fci, £13) plane, where the solid line represents the FP relation as given by equation 
(4) with its observed l-cr dispersion (dashed lines). Bars show the amount of projection effects. Bottom panel: the merging end-products 
in the (fci, £2) plane, where the dashed lines define the region populated by real galaxies as given by BBF. Each model is represented by 
a set of points corresponding to several random projections. Symbols are the same as in Fig. 2. 
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log M../M, t0 

Figure 5. Top panel: stellar central velocity dispersion (normalized to that of the first progenitor) vs. total stellar mass. Points correspond 
to angle-averaged values, bars indicate the range spanned by projection effects. The solid and the two dashed lines represent the Fabcr- 
Jackson relation Lb oc and its scatter, while the dotted line represents Lb °c °o' 5 " Bottom panel: stellar effective radius (normalized 
to that of the first progenitor) vs. total stellar mass. Points and bars have the same meaning as in top panel. The solid line represents 
our "fiducial" Kormendy relation, obtained combining the FP and the FJ with n = 4 (see text); symbols are the same as in Fig. 2. Note 
the different range spanned in the ordinate axes in the two panels. 



scenario (see also NLC02). Models with DM (stars in Fig. 5) 
seem to follow more closely the FJ relation, even though the 
representative points show the same clear trend as the one- 
component merging hierarchy: it seems clear to us that addi- 
tional simulations would bring also two-component models 
completely outside the FJ relation. 

Considering now the agreement with the FP on one 
side, and the failure at reproducing the FJ relation on the 
other side, we would expect that also the Kormendy rela- 
tion is not reproduced by the equal mass merging hierarchy. 
In fact, if we assume that the Es obey to the FJ relation, 
Lb k<7™, then from equation (4) we obtain a fiducial Kor- 
mendy relation (R) e oc L B ' 225_2,/ ". For example, with n — 4 
(for which the exponent is remarkably near to that reported 
by Bernardi et al. 2003) the effective radius of the end- 
product of the last equal mass merger is a factor ~ 3 larger 
than that predicted by the "Kormendy" relation. In other 
words, the results of our simulations are characterized by too 
large effective radii (see Fig. 5, bottom panel). Thus, in the 
explored equal mass merging hierarchies both ao and (R) e 
deviate systematically from the values expected from FJ and 
Kormendy relations. Yet, their deviations from the two scal- 



ing laws curiously compensate reproducing the edge-on FP; 
this shows clearly that the FJ and Kormendy relations are 
not simple projections of the FP, which is actually preserved 
by merging, but, despite their larger scatter, contain addi- 
tional information on galaxy structure and dynamics. 



3.2 Accretion 

In Section 3.1 we showed that, in the case of equal mass 
merging of one and two-component galaxies, the virial ex- 
pectations were qualitatively followed also by the observa- 
tional properties (R) e and ao, in addition to cry and ry. We 
recall that these expectations are independent of the specific 
merging history, as far as the zeroth order seed galaxies are 
identical systems, and so we could consider "accretion" sim- 
ulations superfluous. However, the amount of escaping par- 
ticles or the non-homology of the final end-products could 
depend on the mass ratio of the merging systems, thus re- 
sulting in substantial changes in ay and/or ao. Following 
this argument we now present the results of accretion sim- 
ulations, whose technical setting is described in Section 2. 
From the structural point of view, the end-products of ac- 
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crction are less flattened than the corresponding (i.e. with 
the same mass) equal mass mergers, being, in general, triax- 
ial systems with axis ratios in the range 0.6 < c/a < 0.9, and 
0.7 < b/a < 0.9, and becoming more and more spherical with 
increasing mass as a consequence of the higher number of 
merging events and the random direction of accretion. The 
evolution of the surface brightness profile in the accretion 
hierarchy is shown in Fig. 2, where the best fit Sersic pa- 
rameter m is represented with empty circles and triangles 
(for head-on merging and merging with angular momen- 
tum, respectively). A first relevant difference with the equal 
mass merging results is the decrease of m with mass at mass 
ratios larger than 4 for the head-on accretions, and the con- 
sequent smaller range of m (2 < m < 8) spanned by their 
end-products. Therefore, the explored head-on accretion 
scenario fails at reproducing the relation between the sur- 
face brightness profile shapes and luminosity of real galaxies, 
which are characterized by m increasing with galaxy lumi- 
nosity and also spanning a larger range. A different situation 
is obtained when considering accretion simulations with non 
negligible orbital angular momentum: the structural prop- 
erties of the end-products are curiously similar to the equal 
mass cases. 

In the upper panel of Fig. 3 we plot <tv as a function of 
mass of also of the head-on accretion end-products (empty 
circles). We found that the fraction of mass lost in each 
merging AM/ (Mi + M2), decreases from ~ 0.06 in the first 
step to ~ 0.01 in the last one, summing up to a total amount 
of mass lost over 15 steps of ~ 18%, about the same value 
found in the corresponding equal mass merging cases. The 
slightly larger increase of cry in the head-on accretion hier- 
archy with respect to the equal mass merging cases is fully 
accounted for larger values of v e j , in accordance with the 
estimate given in equation (17). In the lower panel of Fig. 
3 we also show the evolution of tm in the head-on accre- 
tion hierarchy (empty circles). The behavior of this quan- 
tity is qualitatively similar to that of equal mass mergers: 
the slightly flatter slope of the sequence with respect to the 
simple expectation rjj oc M is due to the stronger evolution 
of structural non-homology, as revealed by Fig. 2. On the 
other hand, the end-products of accretion simulations with 
angular momentum (empty triangles in Fig. 3) are more 
similar to equal mass merging case. 

3.2.1 Fundamental Plane 

The most striking difference of the end-products of head- 
on accretion simulations (empty circles) with respect to the 
equal mass merging hierarchy and also to accretion simula- 
tions with angular momentum (empty triangles) is apparent 
in the upper panel of Fig. 4, where the results are plotted 
in the (fci,£;3) space. In fact, after few accretions, the end- 
products are characterized by a £3 decreasing for increasing 
fci, at variance with the FP slope and the trend shown by 
the end-products of equal mass mergers. As a consequence, 
the last explored models (corresponding to an effective mass 
increase of a factor ~ 12) are found at a distance Sks larger 
than the FP scatter. This result, when interpreted by using 
the information in Fig. 2, is not surprising. In fact, being the 
coordinate k.i a measure of non-homology for galaxies with 
constant T* (as those explored in this paper), the decrease 
of the Sersic parameter m at increasing mass (at variance 



with real galaxies) reflects directly in the unrealistic trend in 
the (^1,^3) plane. On the contrary, in the face-on {k\,k2) 
plane (Fig. 4, lower panel), the end-products of accretion 
(both head-on and with angular momentum) evolve along 
the same direction followed by equal mass mergers, albeit 
with a smaller excursion in ki and &2- 

3.2.2 Faber- Jackson and Kormendy relations 

In Fig. 5 we also plot the results of accretion simulations in 
the planes representing the FJ and the Kormendy relations. 
At variance with the case of the edge-on FP, and indepen- 
dently on the presence of angular momentum, the behavior 
of the end-products is considerably similar to that of the 
equal mass mergings, i.e., both FJ and Kormendy relations 
are not reproduced. The end-products have ao smaller and 
{R)c larger than the values predicted by the two scaling re- 
lations for the given luminosity (mass) increase. However, a 
closer inspection of Fig. 5 reveals that both ao and {R) c are 
systematically larger in the case equal mass mergers than 
in the case of head-on accretion: in other words, the end- 
products in the latter case deviate more from the FJ and 
less from the Kormendy than those in the former. The result 
relative to the effective radius simply reflects the different 
evolution of the half-mass radius in the two scenarios (see 
Fig. 3, lower panel), while the behavior of ao depends both 
on structural and dynamical non-homology effects: the com- 
bined effect of the lower values of ao and the larger values of 
{R)e is responsible of the decreasing of ks with ki observed 
in the analysis of the edge-on FP. In other words, at vari- 
ance with equal mass mergers and accretions with angular 
momentum, there is not enough compensation between the 
trends in ao and (R) e able to reproduce the FP tilt in case 
of head-on accretions. 



4 DISSIPATIONLESS MERGING AND THE 
M BH -cr<) RELATION 

On the basis of the results of the previous Section, it is 
particularly interesting to investigate whether dissipation- 
less merging is able to reproduce the Mbh-oo relation, and 
we attempt to answer this question by using the simulations 
presented in Section 3. Here we recall that in our simulations 
we do not take into account the presence BH in the merging 
galaxies, yet we argue that we can reach robust conclusions, 
and this for two reasons. In fact, as well known, at the equi- 
librium (i.e. after the virialization of the end-product) the 
presence of the BH has not significant influence on ao- This 
can easily be seen by considering that the sphere of influence 
of a BH with mass Mbh at the center of a galaxy with central 
velocity dispersion ao has a fiducial radius tbh = GMbh/o"o> 
and combining this equation with equation (7) we obtain 

rB ^°- 01X (2()0I^) Q " kpC ' (W) 

where we assumed Mbh = 1O 8 M when a = 200kms _1 
(see, e.g., Tremaine et al. 2002). The central velocity dis- 
persion used in the definition of the FP and of the Mbh-co 
relations is the luminosity weighted projected velocity dis- 
persion inside the radius {R) c /8 which, for ellipticals with 
ao ^ 200 km s"\ is in the range 0.2 < (R) e /S < 0.6 kpc, 
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Figure 6. Panel a: galactic central velocity dispersion vs. BH mass for classical BH merging; crofi and Mbh.o are the central velocity 
dispersion and BH mass of the first progenitor, respectively. The points correspond to the mean value over the solid angle, while the bars 
indicate the range spanned by projection effects. Solid and dashed lines represent the Mbh-o"o relation for a = 4 and a = 5, respectively, 
while vertical dotted lines show the observed scatter around these best fits. Panel b: same data as in panel a, but for maximally radiative 
BH merging. Panel c: stellar mass vs. BH mass for classical BH merging; M tj o is the stellar mass of the first progenitor and the solid line 
represents the Magorrian (Mbh oc Mb u i gc ) relation. Panel d: same data as panel c, but for maximally radiative BH merging. Symbols 
are the same as in Fig. 2. 



one order of magnitude larger than tbh- Thus, once the nu- 
merical end-products reached the equilibrium, their er is a 
good estimate of the real quantity even in presence of a BH 
of realistic mass. The second reason requires a more careful 
discussion. In fact, during the dynamical evolution of the 
merging system, a BH binary is formed and, in principle, its 
evolution could affect the kinematic properties of the end- 



product even beyond tbh- However Milosavljevic & Merritt 
(2001), with the aid of numerical simulations, showed that 
the formation of a BH binary does not modify significantly 
the central velocity dispersion measured within the standard 
aperture {R) c /8, although the inner density profile is mod- 
ified (for interesting observational evidences of this case see 
Lauer et al. 2002). In any case, we remark here that, if any, 
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the expected effect of binary BHs is to decrease the central 
velocity dispersion, as a consequence of dynamical friction 
heating against background stars. 

On the basis of these considerations we simply assume 
that each seed galaxy contains a BH of mass Mbh,o, and that 
each merging end-product contains a BH obtained by the 
merging of the BHs of the progenitors (at the end of this Sec- 
tion, however, we will discuss two important problem faced 
by this last assumption). Unfortunately, BH merging is still 
a poorly understood physical process, in particular with re- 
spect to the amount of emitted gravitational waves (see, e.g., 
CvA, and references therein), and for this reason we consider 
two extreme situations: the case of classical combination of 
masses (Mbh, 1+2 = Mbh,i + Mbh,2, with no emission as 
gravitational waves) , and the case of maximally efficient ra- 
diative merging (M| H 1+2 = M BH ,i + M BH ,2i corresponding 
to entropy conservation in a merging of two non rotating 
BHs). Following this choice, in Fig. 6 we plot the central 
velocity dispersion of the mergers versus the mass of their 
central BH in the case of classical (panel a) and maximally 
radiative (panel b) BH merging. As expected from the simi- 
larity between the FJ and the Mbh-co relations, in the clas- 
sical case, both (one-component) equal mass mergers (solid 
squares and triangles) and accretion mergers (empty circles 
and triangles) are unable to reproduce the observed relation, 
even when its slope is assumed as high as 5, with representa- 
tive points well outside the observed scatter (vertical dotted 
bars). As it happens for the FJ relation, the situation is 
somewhat better when considering equal mass mergers of 
two-component galaxy models: in presence of DM the stel- 
lar central velocity dispersion is higher for given BH mass 
and the end-products are found closer to the observed rela- 
tion (stars in Fig. 6). Thus, as for the FJ, the reason of the 
failure of dissipationless merging at reproducing the Mbh- 
<to relation is that the end-products are characterized by a 
too low <7o for the given Mbh, i.e., Mbh is too high for the 
resulting cxo- A promising solution to this problem could be 
the emission of some fraction of Mbh as gravitational waves. 
In fact, by assuming maximally efficient radiative BH merg- 
ing, the points are found remarkably closer to the observed 
relation, even if it is apparent that the slope of the Mbh- 
<to relation is not well reproduced by the end-products of 
head-on accretion (empty circles). Note however that, while 
in the classical scenario the Magorrian relation is (obviously) 
nicely reproduced (Fig. 6) , in case of substantial emission of 
gravitational waves the relation between BH mass and bulge 
mass is not reproduced (Fig. 6). 

All the results presented in this Section holds under 
the strong assumption that the BHs of the merging galax- 
ies are retained by the end-products, but there at least two 
basic mechanisms that could be effective in expelling the 
central BHs. The first is related to the general instability 
of three body systems: if a third galaxy is accreted by the 
end-product of a previous merging before the binary BH 
at its center merged in a single BH, then the escape of the 
smallest BH is highly possible. It is clear that if this process 
happens more than a few times, then the Magorrian relation 
will be not preserved at the end. A plausible solution to this 
problem is then to assume that the characteristic time of BH 
merging (first by dynamical friction on the background stars 
and then by emission of gravitational waves) is shorter than 
the characteristic time between two galaxy merging (for a 



detailed discussion of this problem see, e.g., Milosavljevic & 
Merritt 2001, Yu 2002, Haehnelt & Kauffmann 2002, Volon- 
teri, Haardt & Madau 2003). However, a second physical 
mechanism could be more effective in expelling the resulting 
BH from the center of a galaxy merger, i.e., anisotropic emis- 
sion of gravitational waves. This process, commonly known 
as the "kick velocity", is directly related to the fraction of 
BH mass emitted anisotropically during the BH coalescence 
(see, e.g., Flanagan & Hughes 1998, CvA). Gravitational 
waves in fact travel at the velocity of light, and so even 
the anisotropic emission of a few thousandths of the mass of 
the BH binary will produce a recoil (due to linear momen- 
tum conservation) of the resulting BH with a characteristic 
velocity higher than the escape velocity typical of massive 
galaxies. In conclusion, it seems to us not trivial to assume 
that in each galaxy merging the resulting BH will reside at 
the center (see, e.g., Haehnelt & Kauffmann 2000). 



5 DISCUSSION AND CONCLUSIONS 

As outlined in the Introduction, the present study is similar, 
in some aspects, to previous studies concerning the effects of 
merging on the FP. The results of our simulations are consis- 
tent with those of Capelato et al. (1995), based on N-body 
simulations of King (1966) galaxy models, and of Dantas et 
al. (2003), based on Hernquist models, who found that the 
edge-on FP is reproduced by dissipationless merging, ow- 
ing to non-homology of the end-products. However, they 
considered only the first two or three steps of the merg- 
ing hierarchy. In the same line of investigation, Gonzalez- 
Garcia & van Albada (2003) found that also the first gen- 
eration of end-products of dissipationless merging based on 
Jaffe (1983) models preserves the FP when seen edge-on. 
At variance with these works, we used galaxy models with 
a substantially larger number of particles, we explored the 
cumulative effects of several generations of merging, and we 
also compared with other scaling laws, such as the Faber- 
Jackson, the Kormendy, and the Mbh-co relations. 

The main results of our simulations can be summarized 
as follows: 

• The edge-on FP is well reproduced by dissipationless 
hierarchical equal mass merging of one and two-component 
galaxy models, and by accretion simulations with substan- 
tial angular momentum, with their seed galaxies (i.e. the 
merging zeroth order generation) placed on the FP itself. On 
the contrary, in the case of head-on accretion, the k-$ coordi- 
nate decreases at increasing fci, at variance with real galax- 
ies. The physical reason of these different behavior is due 
to a different evolution of structural non-homology of the 
end-products in the two scenarios; in fact the Sersic best fit 
parameter m grows monotonically with mass in equal mass 
mergers (and accretion with angular momentum), while de- 
creases with mass in head-on accretion end-products. 

• In both scenarios the scatter in the (fci, k$) coordinates 
associated to projection effects is of the order of the observed 
1-a dispersion of the edge-on FP. 

• At variance with the results of the edge-on FP, all the 
results of our simulations are in qualitative agreement when 
plotted in the (fci,fc2) plane, where the FP is seen nearly 
face-on. In this plane the mergers move parallel to the line 
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defining the zone of avoidance, and remain in the region of 
the FP populated by real galaxies. 

• As a direct consequence of the numerical findings that 
the "observables" ao and {R) c follow closely the evolution 
of (Jv and ry as predicted by the virial theorem, the end- 
products of our simulations fail to reproduce both the FJ 
and Kormendy relations. However, in the case of equal mass 
merging and accretion merging with angular momentum, 
the large values of (R) c and the nearly constant values of ao 
curiously compensate, and the end-products follow nicely 
the edge-on FP (see the first point above) . 

• Under the reasonable hypothesis that the derived val- 
ues of ao are not strongly affected by the dynamical evo- 
lution of binary BHs, our results show that dissipationless 
merging, while in accordance with the Magorrian relation, 
fails to reproduce the Mbh-0"o relation. We have then shown 
that, allowing for substantial emission of gravitational waves 
during the BHs coalescence, the Mbh-co relation is surpris- 
ingly reproduced, but the Magorrian relation is not. Merging 
by head-on accretion suffers of an additional problem when 
considered together with the Magorrian relation: our sim- 
ulations in fact showed that the best Sersic m parameter 
decreases at increasing mass of the end-product, a behavior 
opposite to what is empirically found (Graham et al. 2001). 

In conclusion the results of the presented simulations 
suggest that substantial dissipationless merging (especially 
when involving stellar systems of very different mass) can- 
not be the basic mechanism of formation of Es. We note, 
however, that our exploration of the parameter space is by 
no means complete, because in the simulations we did not 
consider initial conditions corresponding to multiple merg- 
ing. In any case, indications exist that our results should 
be obtained also under these more general circumstances 
(Nipoti et al. 2003b). Finally, we note that from the numer- 
ical results one cannot exclude, on the basis of FP, FJ, and 
Kormendy relations, that Es could experience few occasional 
mergings with other Es, even in recent times as observed by 
van Dokkum et al. (1999). In any case, it is well known that 
many other astrophysical evidences, based on stellar popula- 
tion properties, such as the Mg 2 -o"o and the color-magnitude 
relations, and the waveband dependence of the relation be- 
tween mass-to-light ratio and galaxy luminosity (see, e.g., 
Pahre et al. 1998b, Pahre, de Carvalho & Djorgovski 1998a), 
strongly argue for a substantial dissipative phase in the for- 
mation of spheroids. 

If we exclude dissipationless merging as the formation 
mechanism for Es, then we are left with two possible solu- 
tions, namely the monolithic scenario or a merging scenario 
where dissipation plays an important role. In this last case 
the presented simulations can also say something. In fact 
gas dissipation should be effective in shrinking the density 
distribution of the merging end-products, and, at the same 
time, increasing their central velocity dispersion, changes 
that are in the direction required by the FJ and Kormendy 
relations. In addition, in order to preserve the Magorrian 
relation, a fraction of the star forming dissipating gas must 
flow on the central BH, producing QSO activity; in this case 
the cosmological evolution of QSOs would be a tracer of the 
merging history, in accordance with the old ages of stars 
in Es. As a consequence, numerical simulations of galaxy 
merging in presence of gas (see, e.g., Bekki 1998), taking 



also into account the feedback from the central BH in a self- 
consistent way (see, e.g., Tabor & Binney 1993, Binney & 
Tabor 1995, Ciotti & Ostriker 1997, 2001), and finally a bet- 
ter understanding of BH merging (see, e.g. Merritt & Ekcrs 
2002; Hughes & Blandford 2002) are highly needed for a 
substantial progress in this field. 
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